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We present an example that shows that the random phase approximation performed on 
high-T^ multi-quasiparticle configurations leads to a rotor picture by calculating excitation 



energies and magnetic properties of 178 W. Then we deduce the effective gn of the high- A" 
, rotors and compare it with that of the low- if one. 

(N 

| Rotation is one of typical collective motions in atomic nuclei. Normally axially 

symmetric nuclei carry large angular momenta in the form of collective rotation about 
ON , an axis perpendicular to the symmetry axis. In some cases in which single particle 

orbitals with large angular momenta jj and their projections to the symmetry axis 
■ Qi sit at the vicinity of the Fermi surface — realized typically in the A ~ 180 region 

— , the nucleus can have large angular momenta by aligning multiple quasiparticles 
(QPs) along the symmetry axis. Sometimes the latter scheme forms yrast states, or 
even if not, isomers are often formed thanks to largeness of K = Qi. 

Detailed information about high-if configurations can be obtained from their 
magnetic properties — static magnetic moments and/or (/-factors inferred from 
B(M1)/B(E2) branching ratios of in-band transitions in rotational bands excited on 
top of high- if configurations. The latter data are transformed into \gx — 9r\/Qo by 
way of the rotor model. ^ Then by assuming appropriate gR and Qo, the extracted 
gx is compared with the weighted average of single- j g-factors with respect to Q%?^ 
On the other hand, (at least lower members of) rotational bands excited on high- 
K configurations can be described as multiple excitations of the precession phonons 
in the language of the random phase approximation (RPA). 3 * 1 ' 4 ) Thus, by calculating 
the wave function of the one phonon state, B{M\ : I = K +1 — > K) can be obtained 
and transformed into the effective (RPA) (gx — gn)', its sign can be determined by 
the calculated E2/M1 mixing ratio. The magnetic moment (//) and accordingly the 

g-factor, g = \ f^(/j) /((J)/^jv)j of high- if configurations can be calculated in the 
mean field level. Since this g essentially coincides with gx, we can deduce gn of the 
high- if rotor by combining the RPA (gx — gn) and g. 

The purpose of the paper is twofold: By applying the above method to 178 W for 
which the richest experimental information 5 ) -7 - 1 is available, we corroborate that the 
RPA gives a rotor picture via the excitation energies and magnetic properties, and 
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deduce gn of high- if rotors. 

We begin with a one-body Hamiltonian, 




Here r = 1 and 2 stand for neutron and proton, respectively, and chemical poten- 
tials A T are determined so as to give correct average particle numbers (N T ). The 
oscillator frequencies are related to the quadrupole deformation parameters £2 and 
7 in the usual way. [We adopt the so-called Lund convention.] The orbital angular 

momentum 1 is defined in the singly-stretched coordinates x' k = \J~^%ki with k = 
1-3 denoting x — z, and the corresponding momenta. Nuclear states with QP 
excitations, i.e., alignments along the x axis, are obtained by exchanging the QP 
energy and wave functions such as 



where ji denotes the signature partner of \i. 

We perform the RPA to the residual pairing plus doubly-stretched quadrupole- 
quadrupole (Q" ■ Q") interaction between QPs. Since we are interested in the preces- 
sion mode that has a definite signature quantum number, a = 1, only two compo- 
nents out of five of the Q"-Q" interaction are relevant. They are K = ±1 components. 
Note that we call the symmetry axis, with respect to which the K quantum number 
is defined, the x axis throughout this paper. That is, we consider the 7 = —120° case. 
These components of the interaction are related to the restoration of the spherical 
symmetry. Requiring the decoupling of this symmetry mode (the Nambu-Goldstone 
mode), J± = J y ± iJ z , the strength of the interaction is determined. Then, utilizing 
the identities in Table III of Ref. 8)*), the RPA equation of motion can be cast into 9 ) 
the following form © , which we use in the actual calculation rather than the original 
equation in terms of Q" operators: 



e^V M ,U M )-»(e^Up,V p ), 



(w 2 - u r 2 ot ) 



A(u) C(u) 
B(u) D(J) 



= 



(3) 



where 



A(u) = u)J y (uj) - LO rot J yz {uj), 
B(uS) = LO rot {J y (uj) - J x ) - uj yz (uj) 

C(u) = Wrot (Jz{oj) ~ Jx) ~ Uj yz (u) 
D(u) = LOj z (uj) - U) rot Jy Z (u)), 



(4) 



*' Strictly speaking, the identities have tiny deviations brought about by the singly-stretched 
(rather than the non-stretched) 1 • s and l 2 potentials. 
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Jz(w) 



Jx = h{J X )/^rot, 
(o=±l/2) 2 

2E^ (iJy(pLV)) 



- E 

(a=±l/2) 

- E 

(a=±l/2) 

E 



(M 2 



2E^ (Jjjw)Y 



E 2 



(5) 



Here we adopt the convention that matrix elements of J y and /iy (below) are pure 
imaginary. The non-spurious part of Eq. (J5J), ^4(w)D(a;) — B{uj)C{uj) = 0, can be 
rewritten as 



ujl eS) (u>)-u Iot [J x -J_ 



r(eff) 



0, (6) 



where the suffixes + and — refer to the AK = +1 and —1 modes, respectively, and 

jf\u)=J 1 _(u)^J yz (u), 



J±{UJ) = Jy 



(7) 



For AK = +1 excitations, corresponding to the precession modes, the excitation 
energy in the laboratory frame is given by 



rot 



rot " 



4 efl V) 



h 2 



(Jx) 



(eflf). 



(8) 



which is independent of a; ro t- Since the excitation energy of the first rotational state 
on the high-ET configuration in the rotor model is given by 



E I=K+1 - E I=K = j {K + 1) 



derived from 



K 2 



(9) 



(10) 



Eq. [{J x ) = K in the cases of 7 = —120° or 60°] and Eq. @ correspond to each 
other well for K 3> 1. In other words, J+\u) in our RPA formalism and J in the 
axially symmetric rotor model correspond to each other. 

Marshalek gave an expression for multipole transition rates, which is valid for 
I ^> 1, in terms of the RPA wave function. 10 ) In the Ml case this reads 

B(Ml:I^I-l) = l([i» y + t i z ,Xl}) 2 , 



11) 
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for the n-th phonon state. Hereafter we concentrate on the precession phonon. By 
equating this with the expression in the rotor model, ^ 

B(M1 :I = K + l^K) = (g K - g R f K 2 (IKW\I - IK) 2 , (12) 

we obtain the RPA \gx — gn\- Its sign is determined by that of the calculated E2/M1 
mixing ratio. 

Calculation is performed for all the high-.K" (4, 6, 8 and 10QP) configurations 
that exhibit rotational bands; K 77 = 13", 14+, 15+, 18", 21", 22", 25+, 28", 29+ , 
30+ and 34+. The K* = 13", 14+ and 15+ are 2v2ir, 18" is 2uAtt, 21" and 22" are 
4^27r, 25 + , 28" and 30 + are 4u4tt, and 29 + and 34+ are 6v4n configurations. 6),7) 
The model space is N osc = 3-7 for neutrons and 2-6 for protons. The strengths of 
the 1 • s and l 2 potentials are taken from Ref. 11). The pairing gaps are assumed to be 
0.5 MeV for 2QP and 0.01 MeV for 4 and 6QP configurations both for neutrons and 
protons. The quadrupole deformation is chosen to be e 2 = 0.235 that reproduces 
in a rough average the value Qq = 7.0 eb that was assumed in the experimental 
analyses. 6 )' 7 ) As for the spin g-factor, gi e{C ^ = 0.7gi iree ^ is adopted as usual. These 
mean that the choice of parameters in this work is semi-quantitative; we checked 
the robustness of the results with respect to their variations. In the cases symmetric 
about the x axis considered here the results do not depend on w ro t, while actual 
calculations are performed at hw m t = 0.001 MeV. 

Figure ^ presents the calculated and observed relative excitation energies of the 
first rotational band members, Ej = k+i — Ej = k- Our RPA calculation reproduces 
their gross features well but with a close look one finds deviations at K n = 18~ , 25 + , 
28", and 29 + that include the vr/ig/2 orbital. Low calculated energies correspond to 
large moments of inertia [see Eq. Q] and their largeness correlates with that of 
calculated Qq. The largeness of Qq indicates the shape polarization effect of this 
high-j orbital to the prolate direction. As for the effect of the vr/ig/ 2 orbital on the 
moment of inertia, see also Refs. 12), 13). 
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Fig. 1. Calculated and experimental excitation energies of the first rotational band members, 
Ei=k-\-i — Ei=k- Data are taken from Refs. 6), 7). 

In FigEJwe compare the RPA (gx — 9r) extracted from the calculated B(M1; I = 
K+l — > K) and the observed one extracted assuming Qq = 7.0 eb from the branching 
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Fig. 2. Calculated intrinsic g of high-JC configurations (solid curve) and calculated (dashed curve) 
and experimental (points with error bars) (<?k — gu). Data are taken from Refs. 6), 7). 



ratios of available lowest transitions in respective rotational bands. Their agreement 
is semi-quantitative. In this figure the calculated intrinsic g = \J~^-{^x) / {{Jx) ^n) of 
high-if configurations are also shown. The two calculated curves roughly correlate. 
According to the relation 1 ) 

K 2 

9R = 9 ~ [9K ~ 9R) jtj + 1 ^ (13) 

with I = K + 1, g almost coincides with gx- Consequently, the difference between 
the two curves essentially corresponds to the effective gR for high-ET cases. Thus, this 
correlation suggests a possibility to deduce the effective gR of the considered high- if 
configurations by substituting the RPA [gx — 9r) t° Eq. Its average value is 

about 0.29 as seen from FigEl This value may be a rough measure of a property 
of high-if rotors. Moreover, an interesting feature is that there are considerable 
variations and those for the configurations including the irhgn orbital are larger 
than others. In order to see it more closely, in Fig|31we compare them with those 
calculated from an approximate relation 14 ) 

9R = rr ~T" 7f 1 ( 14 ) 

Uv \ \J TV 

where the neutron and proton part of the effective inertia, the upper sign of Eq. (JJJ), 
are substituted to J v and J^. It is clear that the gR values deduced from Eqs. lfT5|) 
and (|14[) correspond to each other very well, although those from Eq. (|14|) are much 
larger for the configurations in which the irhg/2 orbital is occupied: The contribution 
to the moment of inertia from the Trhg/2 orbital is large and overestimated in the 
calculation as mentioned before in the case of excitation energies. 

Finally we compare the deduced gR of high-if rotors above and that of the 
ground state band. We calculate g of the 2 + , 4 + , and 6 + states, which to the zeroth 
approximation plays a role of gR for the nearby K ^ configurations, at hw ro t = 
0.053, 0.119, and 0.176 MeV, respectively, with e 2 = 0.235, 7 = 0, and the odd- 
even mass differences A n = 0.883 MeV and A p = 1.026 MeV as pairing gaps. The 
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Fig. 3. gn of high-.fr rotors deduced using Eq. 1131 (solid curve) and calculated using Eq. 1141 
(dashed curve). 

results are g = 0.218, 0.216, and 0.214, respectively, which almost coincides with the 
average for the configurations that do not include the irhgn orbital. This indicates 
that high- IT and low-.fr rotors are similar unless the shape driving irhg/2 orbital is 
included. 

To summarize, we have numerically verified that the random phase approxima- 
tion performed on high-.fr multi-quasiparticle configurations leads to a rotor picture, 
as previously discussed via E2 properties by Andersson et al., 3 ' by calculating ex- 
citation energies and Ml properties. Next we have deduced the effective gR of the 
high-X rotors and compared them with those of the \ow-K rotor near the ground 
state. A more detailed investigation is under progress. 
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